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We propose realizing the quantum anomalous Hall effect by proximity coupling graphene to an 
antiferromagnetic insulator that provides both broken time-reversal symmetry and spin-orbit cou¬ 
pling. We illustrate our idea by performing ab initio calculations for graphene adsorbed on the (111) 
surface of BiFeC> 3 . In this case we find that the proximity-induced exchange field in graphene is 70 
meV, and that a topologically nontrivial band gap is opened by Rashba spin-orbit coupling. The 
size of the gap depends on the separation between the graphene and the thin Him substrate, which 
can be tuned experimentally by applying external pressure. 

PACS numbers: 71.70.Ej, 73.43.Cd, 81.05.Uw 


The possibility of Hall conductance quantization in a 
two-dimensional electron systems in the absence of an ex¬ 
ternal magnetic field was first proposed by Haldane based 
on honeycomb lattice jlj toy model studies. Like the con¬ 
ventional quantum Hall effect, this quantum anomalous 
Hall effect (QAHE) can occur only in two-dimensional 
systems that are not time-reversal invariant. Although 
its appearance is dependent on electronic structure de¬ 
tails, the QAHE is a possibility in any strongly spin- 
orbit coupled magnetic two-dimensional electron sys¬ 
tems. The QAHE has so far been predicted to occur 
in mercury-based quantum wells Q, topological insula¬ 
tor thin films , graphene @-0 and silicene [9] based 
systems, kagome lattice 0 systems, and in thin lay¬ 
ers containing heavy element s 1111 . Related effects can 
occur in optical lattices 0, 0. The QAHE has re¬ 
cently been observed experimentally in chromium-doped 
(Bi,Sb) 2 Te 3 Q. 

Graphene is unique among two-dimensional electron 
systems in its ability to maintain electrical isolation at 
high carrier densities, and is in many ways an ideal sys¬ 
tem for applications of the QAHE. However graphene 
is not magnetic in its pristine form and has extremely 
weak intrinsic spin-orbit coupling lU-17|. In order to 
realize the QAHE effect both spin-orbit coupling and 
magnetism should be induced externally, for example by 
depositing 3d transition metal atoms on graphene |18l |. 
However, metal adatoms tend to nucleate into clusters 
on graphene 0,113. As reported in Refs. 0, ferromag¬ 
netism in topological insulators can be induced due to the 
proximity coupling to ferromagnetic insulators. Compar¬ 
ing to the spreading of surface states of a topological insu¬ 
lator into the thickness direction of a few quintuple layers, 
the single atomic-layer graphene gives a much stronger 
hybridization with the substrate, therefore the proximity 


effect in graphene should be generally larger than that in 
topological insulators. In this Letter we explore the pos¬ 
sibility of depositing graphene on a proper magnetic sub¬ 
strate, which can induce both a time-reversal symmetry 
breaking exchange field and Rashba spin-orbit coupling 
due to broken inversion symmetry. Guided by previous 
wor k@ we expect the QAHE to be realized when the 
substrate satisfies the following criteria: (i) it is insulat¬ 
ing so that graphene provides the only transport chan¬ 
nel; (ii) its orbitals hybridize with those of graphene suf¬ 
ficiently to induce sizable exchange and spin-orbit cou¬ 
pling fields; (iii) its dipolar magnetic field is sufficiently 
weak that the quantum anomalous Hall signature is not 
obscured by normal quantum Hall effect. In this Letter 
we demonstrate that antiferromagnetic insulators with 
properly terminated surfaces can satisfy all these cri¬ 
teria. BiFe0 3 is a heavily studied room-temperature 
multiferroic material (see, e.g. Ref. [22j] and references 
therein). By theoretically depositing graphene on the 
(111) surface, we show that the graphene 7 r bands spin- 
split on a BiFe0 3 thin film substrate. After switching 
on spin-orbit coupling in our electronic structure calcu¬ 
lations, we further find that a small band gap 6 opens 
near the graphene K and K’ Dirac points which are lo¬ 
cated inside the bulk band gap A of the BiFe0 3 thin 
film substrate. Berry curvature analysis confirms that 
the resulting band gap is topologically nontrivial, yield¬ 
ing a two-dimensional quantum anomalous Hall phase. 
Finally, we show that <5 can be increased to a value that 
is large enough for high temperature operation by reduc¬ 
ing the separation between graphene and the thin film 
substrate. 

Our first-principles calculations were performed us¬ 
ing the projected-augmented-wave method 2)| as im¬ 
plemented in the Vienna Ab-initio Simulation Package 
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FIG. 1: (Color online) (a) Supercell for graphene on a (111) 
BiFeOs surface. Gray, purple, yellow, and red colors repre¬ 
sent respectively C, Bi, Fe, and O atoms, (b) A top view 
of graphene with one layer of iron atoms. One Fe atom is 
located at the hollow site labeled by H, while the other two 
Fe atoms are located at the top sites labeled as A and B. (c) 
Band structure along high symmetry lines for a 4 x 4 supercell 
of graphene. 


(VASP) 24|. The generalized gradient approximation 
(GGA) [25| exchange-correlation functional was adopted. 
The value of the Hubbard U was chosen to be 3 eV to 
account for strong intra-atomic interactions on the Fe 
atoms and to correct band gap underestimation. The 
lattice constant of graphene was chosen to be a = 2.39A, 
and the kinetic energy cutoff was set to be 550 eV. Dur¬ 
ing structural relaxation, all atoms were allowed to relax 
along any direction, and forces were converged to less 
than 0.01 eV/A. The first Brillouin-zone integration was 
carried out using 3x3x1 Monkhorst-Pack grids. A vac¬ 
uum buffer space over 15A was included to prevent in¬ 
teraction between adjacent slabs. Hydrogen atoms were 
used to passivate the Bi surface on the side of the BiFeOs 
film opposite to the graphene interface. 

BiFeOs has a perovskite structure in which the mag¬ 
netic Fe ions form a simple cubic lattice. Since BiFeOs 
is a G-type antiferromagnet, the Fe moments in each 
(111) plane are aligned and neighboring (111) planes 
have opposite spin polarizations. We therefore expect 
that a (111) BiFeOs surface, which hybridizes relatively 
strongly with a graphene sheet deposited on it, can in¬ 
duce both a strong homogeneous exchange field and ob¬ 
servable spin-orbit coupling. Because the overall order in 
BiFeOs is antiferromagnetic, it will not induce significant 
magnetostatic magnetic fields in graphene. Figure [T| (a) 
illustrates the supercell used in our calculations, which 
consists of a graphene 4x4 unit cell on a 4-layer BiFeOs 
slab with a (111) surface.The lattice mismatch between 
graphene and BiFeOs in this arrangement is less than 
2%. The average distance between graphene and the 
Fe layer after relaxation is found to be d ~ 2.60A. As 
illustrated in Fig. |T] (b) there are three Fe atoms per 


unit cell in the surface layer: one located at the hollow 
site labeled H and the other two are located at the top 
sites labeled A and B. Since the two top sites belong 
to different sublattices of graphene, the surface structure 
helps to suppress induced staggered sublattice potentials, 
which favors valley-Hall-effect gaps in graphene [2(| and 
are detrimental to QAHE. 

Figure QIc) plots the band structure of an isolated 
graphene layer along the high symmetry lines of the 4x4 
supercell Brillouin zone used in our supercell calcula¬ 
tions. When graphene is placed on top of the BiFeOs 
substrate, it is magnetized by the exchange interaction 
between its 7r-orbitals and the 3d orbitals of the ferro¬ 
magnetic Fe layer, which is antiferromagnetic relative 
to the surface of BiFeOs layer. Figure [2] (a) illustrates 
the spin-polarized band structure of graphene on the 
(111) surface of BiFeOs thin film substrate. Figure [2] 

(b) zooms in the spin-polarized 7r-bands near the Dirac 
point at K’ and shows that the induced exchange splitting 
is M ~ 70 meV. Note that there is a small Dirac point 
splitting between bands of the same spin which can be 
attributed to a small non-zero staggered AB sublattice 
potential in the relaxed structure. We expect that the 
staggered sub lattice potential will be highly random in 
realistic systems and therefore do not expect that this 
splitting will play an important role. Another important 
observation from Fig. [2] (a) is that the Dirac points high¬ 
lighted by the green box are located inside the bulk band 
gap of the thin film substrate. 

When spin-orbit coupling is included a small band gap 
opens at the avoided crossing between crossing bands 
with opposite spin-orientations, as illustrated in Figs. [2] 

(c) and (d). Figure[2](d) zooms in the bands of Fig. Etc) 
near the Dirac point K’ and shows that the avoided cross¬ 
ing gap 5 « 1.1 meV. Experimentally the Fermi level 
Ep can be tuned to lie in this gap by a back gate un¬ 
der the substrate. To determine whether or not the 
resulting insulating state exhibits the QAHE we inte¬ 
grate momentum-space Berry curvatures over the super¬ 
cell Brillouin-zone. Figure [3] plots the Berry curvature 
f l(k) summed over occupied valence bands: 


Q{k) 


n n'^n 
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(1) 

along high symmetry lines. Here v x and v y are operator 
components along the x and y directions and f n = 1 for 
occupied bands. One can see that the Berry curvature is 
non-zero and that it has the same sign near the K and K’ 
valley points. As a consequence its integration over the 
Brillouin zone must give rise to a nonzero Chern number 


C: 



( 2 ) 


The quantized anomalous Hall conductance a xy is related 










3 



FIG. 2: (Color online) Bulk band structure of graphene on BiFeOs. (a) Spin-resolved band structures, (b) Zoom-in on the 
green box of panel (a). The blue and red in panels (a) and (b) plot the spin-up and spin-down bands, respectively, (c) Band 
structure when spin-orbit coupling included. A bulk gap <5 of about 1.1 meV is centered on the avoided crossings near the K 
and K’ valleys. Note that the nontrivial band gap 5 is inside the bulk gap A of BiFeOa. (d) Zoom in on the green box of panel 
(c). 



FIG. 3: (Color online) Berry curvature along high symmetry 
lines summed over occupied valence bands below. 


to the Chern number C by a xy = Ce 2 /h. Due to the ab¬ 
sence of external magnetic field, our obtained insulating 
state is the long-sought quantum anomalous Hall effect. 

According to Ref. [6j, both engineered exchange¬ 
coupling and engineered Rashba spin-orbit coupling must 
play a crucial role in any graphene based QAHE. In the 
present case, both terms arise from the interaction be¬ 
tween graphene and the BiFeC >3 thin film substrate. The 
former originates from the proximity exchange coupling 
and the latter by breaking the mirror symmetry of the 
graphene plane. A staggered sublattice potential can also 
influence the quantum anomalous Hall gap. The low- 
energy continuum model of graphene when all of these 
effects are present is: [$, [Tfl 


H(k) 


vf (r/cp,; k x T T A/T aSz 

T &y&x') T U(T z l s , 


( 3 ) 


where Vf = 3t/2 is the Fermi velocity, t ~ 2.60 eV, 
r) = ±1 for K , K ’ respectively, and er and s are Pauli 
matrices that act on sublattice and spin degrees of free¬ 


dom. The first term in Eq. [3] is the Hamiltonian of pris¬ 
tine graphene, the second describes 7r-electron exchange 
coupling to the magnetic order, the third is Rashba spin- 
orbit coupling and the final term is the staggered sub lat¬ 
tice potential. By fitting Eq.[3]to the first-principles band 
structure near the Dirac points, we can extract values 
of the three coupling parameters. For example, Rashba 
spin-orbit coupling in the present system is estimated to 
have coupling constant A^j ss 1.26 meV. 

Since both exchange field and Rashba spin-orbit cou¬ 
pling are induced by hybridization between graphene and 
thin film substrate orbitals, we expect that they will 
strengthen if the separation between the graphene sheet 
and the thin film substrate is reduced. In experiment 
this can be achieved by applying an external uniaxial 
stress perpendicular to the graphene sheet. In Fig. [4] 
we plot the band gap 5, and the exchange field M and 
Rashba spin-orbit coupling Xr strengths vs. separation 
change do — d. We find that Xr increases by a factor of 
~ 3, M by a factor of ~ 2, and 5 by a factor of ~ 3 
when the separation is reduced by ~ 20%. As expected, 
the band gap 6 increases with decreasing d [square sym¬ 
bols in Fig. [Ha)]. When d = 2.15 A, the band gap S 
can reach 3.50 meV, corresponding to a temperature of 
40.6 K. The low-energy model band gaps in Figure (Ha) 
(red circles) do not account for any potential influence of 
the staggered sublattice potential. The close agreement 
demonstrates that this part of the Hamiltonian does not 
play a major role in determining the size of the gap. 

In Figure [5] we plot the intrinsic [2!|| anomalous Hall 
conductivity. These results were calculated by integrat¬ 
ing Berry curvatures calculated using a lattice represen¬ 
tation of Eq. 0 over the Brillouin-zone for M = 142 
meV and Xr = 4.2 meV, corresponding to (do — d) = 
0.45 A. As expected the Hall conductivity is quantized 
(a xy = —2e 2 /h ) when the Fermi level lies inside the bulk 
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FIG. 4: (Color online) (a) Bulk band gap 5 as a function 
of a change (do — d) in the separation between the graphene 
sheet and the BiFeOs (do is the optimal distance obtained 
from the full relaxation). The square and circle symbols refer 
respectively to gaps from density ab initio calculations and 
gaps from the low-energy model with the staggered sublattice 
potential set to zero, (b) Exchange field M as a function of 
(do — d). (c) Rashba spin-orbit coupling strength as a function 
of (do — d). 



FIG. 5: (Color online) Anomalous Hall conductance as a func¬ 
tion of Fermi level calculated from the tight-binding model 
Hamiltonian corresponding to Eq. (El- 


banc! gap. We find that the intrinsic anomalous Hall con¬ 
ductivity rapidly approaches zero when the Fermi level 
is outside the gap. This is expected since the weak 
Rashba spin-orbit coupling implies that only bands near 
the band gap contribute to the nonzero Berry curva¬ 
ture [see Fig. [3]. 

In summary, by performing explicit calculations for 
graphene on the (111) surface of BiFeC >3 we have demon¬ 
strated that the QAHE can be realized in graphene. We 


expect that our results apply to graphene on any un¬ 
compensated antiferromagnetic insulator surface. Since 
many insulating materials have antiferromagnetic order 
at room temperature, this effect raises the prospects for 
a high temperature QAHE which might have interest¬ 
ing potential for applications. In the BiFeCU case the 
topologically non-trivial gap is not sufficiently large to 
achieve this goal, although it can be enhanced by ap¬ 
plying pressure. Our work motivates a search[30j for 
other substrates. For example, the (100) surfaces of 
A-type antiferromagnets are uncompensated and would 
work if the surface moments have a nonzero component 
along the surface normal. When the van de Waals’ 
interaction is further considered, the distance between 
graphene and the semi-infinite substrate is slightly de¬ 
creased to do = 2.58A, which naturally increases the re¬ 
sulting Rashba spin-orbit coupling and the QAHE gap. 
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